Impact of mergers on LISA parameter estimation for nonspinning black hole binaries 
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We investigate the precision with which the parameters describing the characteristics and location 
of nonspinning black hole binaries can be measured with the Laser Interferometer Space Antenna 
(LISA). By using complete waveforms including the inspiral, merger and ringdown portions of the 
signals, we find that LISA will have far greater precision than previous estimates for nonspinning 
mergers that ignored the merger and ringdown. Our analysis covers nonspinning waveforms with 
moderate mass ratios, q> 1/10, and total masses 10 5 < M/M@ < 10 7 . We compare the parameter 
uncertainties using the Fisher matrix formalism, and establish the significance of mass asymmetry 
and higher-order content to the predicted parameter uncertainties resulting from inclusion of the 
merger. In real-time observations, the later parts of the signal lead to significant improvements in 
sky-position precision in the last hours and even the final minutes of observation. For comparable- 
mass systems with total mass M/Mq ~ 10 6 , we find that the increased precision resulting from 
including the merger is comparable to the increase in signal-to-noise ratio. For the most precise 
systems under investigation, half can be localized to within O(10 arcmin), and 10% can be localized 
to within 0(1 arcmin). 
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I. INTRODUCTION 

Gravitational waves carry a tremendous amount of in- 
formation through the universe. It is the goal of the 
emerging field of gravitational wave astronomy to access 
that information and bring it to bear on the problems of 
astrophysics and cosmology. The current generation of 
gravitational wave detectors, such as the Laser Interfer- 
ometer Gravitational-wave Observatory (LIGO) 0], are 
focused on the detection of gravitational waves from iso- 
lated astrophysical systems. LIGO and its contempo- 
raries will also provide some minimal information on the 
parameters of these systems such as their mass, luminos- 
ity distance, and approximate sky position. The quality 
of this information will be limited by the relatively low 
signal-to-noise ratios (SNR) expected for LIGO events 



The Laser Interferometer Space Antenna (LISA) || , a 
space-based detector of gravitational waves in the milli- 
Hertz band, will detect the inspiral and merger of super- 
massive black hole binaries (BHBs) with very large SNRs 
(100 - 10 4 ) out to redshifts of z ~ 10 or greater 0,0. 
These large SNRs make it possible to extract a large 
amount of information from each event including mass, 
mass ratio, spins, orientation, luminosity distance, and 
sky position. Because sources of gravitational waves are 
strongly dominated by gravitational dynamics, and be- 
cause the waves are expected to propagate through inter- 
vening matter with little interaction, these observations 
may provide an unusually clean and direct measurement 
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of the source system parameters. Of particular inter- 
est are the distance and sky position parameters, which 
will drive LISA's ability to narrow the set of candidate 
source galaxies or clusters for merger events, potentially 
opening up a range of multi-messenger astronomy op- 
portunities. For instance, the coincident measurement of 
gravitational and electromagnetic signatures from a sin- 
gle source galaxy (i. e. standard sirens 0) will allow a di- 
rect measurement of the redshift-luminosity relationship, 
thereby constraining the dark-energy equation of state. 
While cosmological models predict that the dark-energy- 
dominated era began fairly recently at z ~ 1, measure- 
ments for much larger redshifts are not currently possi- 
ble by other methods, and the standard siren method is 
limited only by the achievable range of coincident elec- 
tromagnetic observation. 

Extracting this information requires a waveform model 
that can provide templates with sufficient fidelity to dis- 
tinguish between signals with different parameters in the 
presence of instrumental noise. For BHBs, the complete 
waveform signal is traditionally divided into three differ- 
ent regimes: the inspiral, which can be described using 
post-Newtonian (PN) orbital dynamics; the ring-down, 
which can be treated using black-hole perturbation the- 
ory; and the merger, which bridges the two and can be 
predicted using numerical relativity. 

Ideally, an estimate of LISA's ability to measure the 
parameters of observed BHBs would include the informa- 
tion contained in the complete waveform. However, the 
difficulty associated with modeling the merger has led 
the majority of such studies (the exceptions being Q, 
3, and @) to include only the inspiral portion of the 
waveform. Until recent advances in numerical relativity 
[l0l - [l2j opened the door to a complete understanding of 
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General Relativity's predictions for these signals, it was 
not clear whether theoretical knowledge about the final 
strongest moments of the signal would be available for 
system parameter measurements. A naive guess at the 
consequences of omitting the merger would be that the 
loss in parameter precision would be proportional to the 
loss in SNR. This assumes that the two portions of the 
signal have equal density of information per unit SNR 
and that that information is independent. 

There are two reasons that are sometimes cited for ex- 
pecting that the effect of the merger on parameter preci- 
sion will be less than that on SNR. The first is that the 
merger encompasses very few GW cycles compared with 
the observed portion of the inspiral and it is expected 
that information content correlates with the number of 
cycles. The second is that, in the low-frequency limit, 
the sensitivity of LISA to parameters such as sky posi- 
tion is entirely generated by the orbital motion of the 
LISA constellation and the merger is too short in dura- 
tion to experience a significant orbital modulation. 

In this work we investigate the significance of the 
merger to LISA parameter estimation using quasi- 
analytic waveforms that are tuned to match the results 
of numerical relativity. We restrict ourselves to non- 
spinning binaries with moderate mass ratios (g<l/10) 
and explore the parameter space around a candidate sys- 
tem with total redshifted mass of 1.33 x 1O 6 M0, mass 
ratio q = 1/2, and redshift z = 1. Astrophysically, we 
expect black holes to have spin. While including spin 
fl3l fli| and mergers each separately improve parameter 
estimation, it is not known how these effects will com- 
bine. Therefore, including both spin effects and mergers 
will be an important followup to this investigation. 

In section|nl we discuss the methods employed for gen- 
crating models of the complete waveform signals and the 
instrument, and for estimating parameter measurement 
precision. In section IIII1 we examine LISA's ability to 
measure binary black-hole system parameters. The pri- 
mary novelty of our results is that we assume theoretical 
knowledge of the complete signal is applied in the obser- 
vational analysis. We examine the impact of including 
the merger and higher harmonic content (jlll A[) for com- 
parable mass systems near 10 6 M Q , the variation of the 
results across a range of masses ([III B[) and mass ratios 
([III C[) . In pil E[) we study how sky position information 
accumulates in time, which will impact how LISA ul- 
timately interacts with other astronomical instruments. 
We summarize our key results in section IIVI 



II. METHODOLOGY 

Before presenting our findings, we will briefly review 
the steps taken to estimate the precision with which LISA 
will measure astrophysical parameters using complete 
waveforms. Theoretical predictions of the strong-field 
gravitational dynamics and radiation generation must be 
encoded in a parameterized waveform signal model. We 



then need to apply a model of the instrument, including 
the response to signals and the sources of noise. Finally, 
we need to estimate the theoretical limit on the uncer- 
tainty of the measured signal parameters that could be 
achieved from a measurement consisting of our realiza- 
tions of the signal and noise content. 



A. Waveform model 

We assume that Einstein's theory of gravity correctly 
describes black hole binary systems. While numerical rel- 
ativity can now treat the final moments of these events, 
it would be impractical to conduct simulations covering 
the parameter space of interest. Furthermore, general 
signal templates must cover the complete signal includ- 
ing the long-lasting inspiral signal which can not be mod- 
eled numerically, but which is well-described by the post- 
Newtonian (PN) approximation (l5j . Instead our com- 
plete waveform model is based on a variation of the PN 
treatment that has been "tuned" to approximate the nu- 
merical simulation results at late times. Such a model can 
be tuned using available numerical data, while providing 
reasonable, if unverified, model signals for arbitrary pa- 
rameter sets. 

In order to investigate LISA's capabilities for recov- 
ering source parameters, we specifically use a waveform 
model [l6[ tuned to match the available numerical simu- 
lations for nonspinning black hole binaries. This model, 
referred to as the IRS-EOB model, uses a conventional 
effective-one-body (EOB) Hamiltonian formalism for the 
adiabatic inspiral [17| . For the merger-ringdown, a fit to 
a physically-motivated functional form is employed for 
the phasing (see Eq. 9 in [IH), while the amplitude is 
calculated using a model for the flux that is constrained 
both to be consistent with the inspiral flux through 3.5 
PN order, and also to vanish as it approaches the ring- 
down frequency (referred to as "Model 2" and given by 
Eq. 19 in [IH). The physical motivation, that the radia- 
tion can be treated as though it were being generated by 
a shrinking rigid rotator, explains the IRS in IRS-EOB, 
which stands for "implicit rotating source" . 

For a unit mass-system (m% + to 2 = 1), the source 
model depends only on the remaining intrinsic source pa- 
rameters, the mass ratio q = mxjm^ (where mi > mi), 
and the spins, which we set to vanish for this initial 
investigation. Throughout this work, we employ wave- 
forms that correspond to ~ 10 6 M of observation, or ~ 3 
months for our fiducial case with M = 1.33 x 10 6 M Q . 
Here, we use units of G = c = 1, so that 1M = 
4.92 x 10~ 6 (M/Mq) seconds. Therefore, lower -mass sys- 
tems require longer simulations in M. We are limited 
computationally from employing longer waveforms, but 
we have verified that our results do not change signif- 
icantly (with the single exception of the uncertainty in 
total system mass) by doubling the waveform length for 
the lowest-mass cases investigated. We use a model ca- 
dence of 0.5 M. The signal is resampled when we apply 
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the detector's response function, so that the final signal 
cadence corresponds to a quarter wavelength at the high- 
est frequency reached by the I = 4, m = +4 harmonics. 
This is true even for cases where we restrict the calcula- 
tion to have only quadrupolar content. 

After the source calculation, we derive the incident 
waveforms referenced to the solar system barycenter 
(SSB), at which point we can apply the response of the 
LISA detector. Computation of the incident waveform 
in SSB frame depends on eight additional parameters: 
the redshifted total system mass M — M a (l + z) (with 
M Q the rest mass and z the redshift), luminosity distance 
Dl, coalescence time t c , and three angles describing the 
orientation of the binary, for which we use the inclina- 
tion i (using the convention that i = corresponds to 
the line of sight being coincident with the orbital axis 
of the binary), initial orbital phase <j) a and the polariza- 
tion phase ip At the source, the emitted radiation can be 
decomposed in spin-weighted spherical harmonic compo- 
nents hi m of the dimcnsionless gravitational wave strain 
(scaled for unit distance from the source). Here the 
strain is complex-valued to represent both polarization 
components, h = h+ + ih x . Specifying the parameters 
(t, 0Oi M, Dl) allows us to calculate the solar-system 
incident waveform hs'- 
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where ~~ Yg m are the spin-weight —2 spherical harmonics 
[lsj . The two additional parameters, the ecliptic latitude 
P and longitude A, describe the sky location of the binary 
in the SSB frame. The dependence on sky location is ap- 
plied by the instrument response, which we discuss below. 
We use the vector A a = (InM, InD^, f3, A, i, (f> , ip, t c ) to 
denote the complete set of variable parameters. Note 
that the dependence on mass ratio, q, is not explicitly 
included in A a but is instead implicitly included in the 
hi m . This is because q is not varied when computing 
parameter uncertainties (see section III CP , a procedure 
consistent with that used in Q . This is equivalent to the 
assumption that there is no uncertainty in the measure- 
ment of q. We intend to relax this assumption in future 
investigations. 



B. Instrument model 

The instrument model consists of two components: a 
prescription for converting hs into signals observed by 
the instrument, and a description of the instrument noise. 



The LISA instrument consists of a constellation of three 
spacecraft located at the vertices of an approximately 
equilateral triangle with a side length of 5 x 10 9 m. The 
light travel time along each of the six one-way links is 
monitored using laser interferometry. These individual 
link measurements are then combined using a technique 
known as Time Delay Interferometry (TDI) [l9[ to yield 
observables that contain gravitational wave signals and 
suppress instrumental noise. Of the many families of TDI 
observables [2(|, the ones most suitable for data anal- 
ysis are the orthogonalized or "optimal" variables plj . 
For this investigation, we have developed a set of orthog- 
onal variables we refer to as "pseudo-A, E, T" (here- 
after A = {A, E,T}), which are analogous to the origi- 
nal A = {A, E, T} variables in [2l[ except that they are 
constructed from the Michelson X = {X, Y, Z} variables 
rather than the TDI generators a = {a, f3, 7}. The soft- 
ware package Synthetic LISA [22j is used to generate the 
X variables from the incident gravitational waveforms, 
hs and the A variables are then computed as 
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(2) 



An overall factor of | has been applied to all three for- 
mulas to make A and E agree with A and E in the low- 
frequency limit. 

Synthetic LISA can also be used to model instrument 
noise. However, Synthetic LISA includes statistical fluc- 
tuations in its noise-generation algorithms, whereas we 
are interested in studying parameter uncertainties that 
result from differences in the incident waveforms for typ- 
ical instrumental noise levels. We therefore wish to sup- 
press the impact of statistical fluctuations in a given re- 
alization of the noise. One way to do this is to generate 
an ensemble of noise realizations and average them. Sat- 
isfactory results can be achieved with a suitable number 
of averages at the expense of increased computational 
effort. For this analysis, we have followed the proce- 
dure used in [23[ to produce analytic estimates of the 
mean power spectral densities of the noise in the TDI 
channels directly from the acceleration and optical path 
length noises in the individual links. This procedure re- 
quires making assumptions such as stationary, equal arm 
lengths. The expressions for the one-sided spectral den- 
sities in the A observables arc 
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S A ,B = 2sin 2 ($)[2(3 + 2cos($)+cos(2$))S'p m + (2 + cos($))S , p] 
S T - 8sin 2 ($)sin 2 ($/2) [4sin 2 ($/2)S pra + S op ] , 



(3) 
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where $ = loL/c and L is the arm length, expressed as 
a light-travel time. These expressions are the analog of 
Eqs. 67 and 68 in and Eqs. 19 and 20 in [H| for the 
noise response of the original A, and we have verified that 
we duplicate the results in [2l[ for A using a (accounting 
for a typographical error which appears in Eq. 20 of 21 1 
which, if corrected, would make it consistent with 24 1 
and with our results). 

The quantities S pm and S op are the one-sided spectral 
densities of the proof mass acceleration and optical path- 
length noises, respectively, expressed as equivalent strain. 
They are modeled as 
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The acceleration noise power spectrum in S pm includes 
an additional / _1 reddening below 0.1 mHz to account 
for the unmodeled behavior of the instrument below the 
LISA band. 

As a check of the expressions in Eq. [3J we used Syn- 
thetic LISA to model an ensemble of 1000 realizations of 
noise in the TDI A channels. Fig. Q] shows a compari- 
son of the mean power spectra of this ensemble with the 
expressions in Eq. [3J In general, the agreement between 
the simulated noise and the analytic expressions is quite 
good. Deviations between the two curves indicate areas 
of potential concern when evaluating SNR and parameter 
sensitivity. For example, the analytic noise expressions 
in Eq. [3] contain nulls at frequencies corresponding to the 
inverse round-trip times of the constellation. The simu- 
lated data, be it signal or noise, is finite at these frequen- 
cies due to spectral estimation effects. This leads to a 
spurious divergent contribution to the SNR and param- 
eter uncertainties at these frequencies. To guard against 
this, we applied a noise floor of 10 _40 (//1 Hz) 2 , elimi- 
nating the nulls in the noise response. In addition the 
'flexing' of the LISA arms due to orbital variations was 
disabled in Synthetic LISA to maintain consistency with 
the expressions in Eq. [3J which were derived assuming 
constant arm lengths. 



FIG. 1: Comparison of different noise realizations for the A, 
E, T channels. We have employed a noise model (light dash- 
dotted for A and E, light dotted for T), and verified that 
it agrees well with an averaged ensemble of simulated noise 
using Synthetic LISA (dark solid for A and E, dark dashed for 
T). Galactic foreground noise is not included in these traces 
but is included in parameter sensitivity calculations. 



in the remainder of our analysis. We note that, since the 
information from T will be present at high frequencies, 
its exclusion will lead us to produce conservative uncer- 
tainty estimates with systematically worse uncertainties 
than might be otherwise obtained for higher-mass cases 
where T-channel response can be non-negligible. Our 
treatment is consistent with previous studies which have 
generally neglected the details of LISA high-frequency 
response. We plan to include T in future work. 

The final component of the noise model is the fore- 
ground of gravitational waves from unresolved compact 
binaries. We use the model for the galactic foreground 
that was developed in (25|. Specifically we add to the 
expressions for g in Eq.[3Ja galactic foreground noise, 
Sgai, given by 



5 gaI = [4$sin($)] 2 S, 



conf 5 



(5) 



where S^onf is taken from Eq. 14 of [25| . The contribution 
from S'gai is not included in Fig. Q] but is included in all 
SNR and parameter estimation calculations. 



The other of disagreement between the simulated 
and analytic noise in Fig. Q] is at the low-frequency end 
of the T channel. The analytic expression in Eq. [3J pre- 
dicts that the noise in the T channel should continue 
to decrease with decreasing frequency, while the simu- 
lated noise levels off. As a result of using Synthetic LISA 
to model the signal response, and Eq. [3J to model the 
noise, there would be a large spurious contribution to 
SNR and parameter sensitivity at low frequencies, pre- 
cisely the band where the T channel is not expected to 
contribute. As the source of this discrepancy has not yet 
been identified, we have elected to exclude the T channel 



C. Parameter estimation using the Fisher matrix 

To approximate the measurement precision that LISA 
can achieve, one approach we can take is the Fisher ma- 
trix formalism. If the LISA data stream consists of a 
waveform, h(A a ), embedded in a signal, s, so that the 
noise, n, is given by n = s — h, then the probability that 
a signal contains a waveform with the parameter set A a 
is given by the likelihood function, 

p(A a \s) OC e -(HA")s\h(A«)-s)/2 (6) 
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where (• • • | ■ • • ) is a noise- weighted inner product [26| . 
The "maximum likelihood" set of parameters, A a , is the 
one that maximizes p. Errors in the A a set of parameters 
can be assessed by expanding p around A a , such that 

p(A a \s) cx e -r^A^AV2 (7) 

where 8A a = A a — A a . The Fisher information matrix, 
r a 6, which is the centerpiece of our subsequent analysis, 
is defined to be 

where a and b are parameter indices. Throughout this 
work, we calculate the parameter derivatives in Eq. [5] 
using one-sided differencing, with a fractional step size 
e a = AA a , where we set A = 10~ 4 for the coalescence 
time, and A = 10~ 6 for all other parameters. 

To lowest order in an expansion in SNR -1 , the covari- 
ance matrix, E ab , is just the inverse of the Fisher matrix: 

E ab = (pafc^ 1 ^ + ^(SNR- 1 )] , (9) 

so that <j a = VS aa is the standard deviation of param- 
eter a. The covariance matrix is symmetric, with the 
off-diagonal terms giving the covariance between param- 
eters, and the diagonal terms giving the variance of each 
parameter. Because inverting the Fisher matrix to find 
the covariance matrix is not always valid, we verify our 
results by testing individual cases at random for each 
system of interest using Markov Chain Monte Carlo sim- 
ulations. 

Because we use InM and ln£>£ as parameters, the re- 
sulting uncertainties are fractional: 

a lnM w a Mj Mt 
a lnD L ^ a D L/DL _ (10) 

We therefore express these as fractional uncertainties 
throughout this work. Though we will refer to the quan- 
tities a a as "uncertainties" throughout this work, we wish 
to note that they are a measure of precision, not neces- 
sarily accuracy. 

Another parameter of interest is the precision of the 
sky localization, expressed as the area of the uncertainty 
ellipse on the sky, Q. This is sometimes referred to as 
Afl or cr n in the literature, but we simply use f2, given 
that it is a measure of an area of uncertainty, rather than 
a measure of uncertainty of an area. To construct this 
value from our data, we make the approximation 

ft = 2tt\J((j x a c ° s P f - (£A : co S/ 3) 2 (ii) 

« 27rsin/3y / (cr A af 3 f - (E A ^) 2 . 

We note that generally SI ^ ir a x a 13 , both because we 
quote median values throughout, and because we define 
the ellipse as in [27| , so that the effect of Eq. [11] is to 
diagonalize the Fisher sub-matrix formed by /3 and A, 
and thereby calculate the area using the true semi-major 
and semi-minor axes. 



D. Information Accumulation 

One aspect of LISA's parameter sensitivity that is 
of interest is the way in which parameter uncertainty 
evolves with time. A simple way to investigate this is 
to truncate the signal at a specified sequence of times 
before merger. To avoid edge effects in the ensuing spec- 
tral estimation, the truncated signal is tapered with a 
raised cosine window with a length of approximately 
one wave cycle of the quadrupole mode at the time of 
truncation. The Fisher matrix and covariance matrix 
are then computed using this truncated signal. By con- 
structing a sequence with progressively later truncation 
times, one can trace the evolution of parameter uncer- 
tainty. Fig.[2]shows the time-evolution of the uncertainty 
in the ecliptic latitude, /3, for an equal-mass system with 
M = 1.33 x 10 6 M o at z = 1. Three variants of the 
time curves are shown corresponding to different taper- 
ing conventions. The taper either starts at the designated 
truncation time, ends at that time, or the mid-point of 
the taper occurs at that time. This distinction, while 
seemingly trivial, shows the potential impact of slight 
differences in the treatment of time-domain signals with 
regard to taper length, type, or placement, as clear dif- 
ferences can be seen in Fig. [21 For early times, the taper 
prescription does not matter. At later times, the precise 
taper implementation becomes more important, leading 
to a time shift between the three curves. In our discussion 
of information accumulation in section IIIIE1 truncation 
time refers to the mid-point of the taper. 

For signals comprising a single harmonic (or rather, 
a single pair of £, +to modes, which are complex con- 
jugates of each other for nonspinning waveforms), it is 
straightforward to repeat this analysis in the frequency 
domain. This provides an internal consistency check for 
our code. For each desired truncation time, the instanta- 
neous frequency of the source waveform is used to com- 
pute the corresponding signal frequency. The evaluation 
of the inner product in Eq.[8]is then limited to frequencies 
below this. The curve in Fig. [5] shows that this approach 
is consistent with the time-domain approach, at least up 
to times very near the merger. It would be possible to ex- 
tend this technique by tracking each mode separately and 
computing a different frequency cutoff for each. However, 
we find the time-domain approach to be more straight- 
forward. 



E. Caveats 

It is well known that the Fisher-matrix approach is 
prone to a number of potential pitfalls (28[. In this sec- 
tion, we attempt to address a few of them. 

The first potential issue is our approximation of the 
parameter derivatives, dh/dA a , using a one-sided finite 
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FIG. 2: Estimated improvement in ecliptic latitude uncer- 
tainty as a function of signal truncation time, t, relative 
to the coalescence time, t c for an equal-mass merger with 
M = 1.33 x 10 6 M Q at z = 1. The curves labeled "pre-taper", 
"mid-taper", and "post-taper" were computed using time- 
domain truncation, with the truncation time corresponding 
to the end of, middle of, and beginning of a one-radiation- 
cycle-long taper, respectively. The curve labeled "frequency" 
was computed using the full waveform but imposing an upper 
frequency cutoff when evaluating the inner product in Eq. [8] 
corresponding to the instantaneous signal frequencies at the 
truncation times. The vertical dash-dotted line corresponds 
to the Schwarzschild ISCO (see Sec. lIII A|) . The vertical thick 
solid line separates the times prior to merger from the times 
after merger. 



M 


Dt 


R 
H 


\ 




ib 








4.8e-2 


5.7e-2 


1.3e-2 


5.0e-2 


4 Qe-2 4 Se-2 


6.8e- 


4 


4.4e- 


2 




3 0e-5 


1.9e-4 


2.6e-5 


1 np-S Q 1p-4 


1.9e- 


4 


4 0e- 


5 






1.4e-3 


1.3e-3 


4.0e-4 2.1e-4 


3.9e- 


4 


1.5e- 


4 








l.oe-4 


o.ze-o /.le-o 


l.oe- 


4 


o.oe- 


- 












5.1e-5 l.le-2 


1.3e- 


3 


1.3e- 


4 










1.6e-4 


2.1e- 


4 


7.1e- 


3 












4.0e- 


4 


5.2e- 


5 
















3.4e- 


4 



TABLE I: Consistency of the covariance matrix for an equal- 
mass system with M = 1.33 x 10 6 M at z = 1. The tabu- 
lated quantity is the symmetric matrix |1 — E £1 /E E2 |, where 
E Ei is a covariance matrix calculated using the Fisher matrix 
method with a one-sided finite-difference step size of For 
this comparison, si = 10 -6 and £2 = 10 -7 . 



signals throughout the relevant parameter space is likely 
to be available at the time of LISA's operation. We there- 
fore focus on the additional source information that may 
be obtainable when the full signal predictions are applied 
in the observational analysis. 



III. RESULTS 



difference approach. As a rough check of the validity of 
this approximation, we have computed the Fisher and 
covariance matrices using various finite difference step 
sizes and verified that the results were consistent. Table 
H] demonstrates that an order of magnitude decrease in 
the step size changes covariance terms by a few percent 
at worst. 

A related concern is that the Fisher information ma- 
trix precision estimates assume that the relevant portion 
of the likelihood function can be treated as a quadratic 
function. This assumption should be guaranteed by the 
large SNR, but is not explicitly verified. 

Finally, actual observations will require the implemen- 
tation of concrete algorithms for exploring the likelihood 
function, such as those pursued in the Mock LISA Data 
Challenges [29]. In those Challenges, it has been demon- 
strated that accuracy may be impacted, for instance, by 
complicated structure in the likelihood function, includ- 
ing multiple maxima and extended shallow regions. Sys- 
tematic errors are also possible, for instance, if errors in 
the theoretical signal predictions should exceed statisti- 
cal errors [301 ] . Highly accurate merger waveform predic- 
tions, and corresponding models tuned to those predic- 
tions, are currently available only for a very limited sam- 
pling of specific black hole system configurations. How- 
ever, this area of study is advancing rapidly, and it now 
appears that accurate information about the complete 



Since there are variations in the parameter uncertain- 
ties across the parameter space, we perform Monte Carlo 
simulations to find the distribution of uncertainties For 
specific choices of masses and luminosity distance (pa- 
rameters M, q, and Dl), we conducted Monte Carlo sim- 
ulations consisting of 1024 randomly-generated parame- 
ter sets for all the cases shown, and have spot-checked 
that our results do not change significantly if we increase 
to 8192 parameter sets. The remaining parameters are 
drawn from uniform distributions, with t drawn from a 
uniform distribution in cos 1 to give uniform sky coverage. 



A. Adding the merger and higher harmonics 

Table [H] summarizes the improvement in parameter 
uncertainties resulting from the addition of merger for 
the cases of equal-mass systems and mass-ratio q = 1/2, 
each with a total mass of 1.33 x 10 6 M Q at a rcdshift 
z = 1. We compare the uncertainty estimates obtained 
with four different options for the waveform models. Two 
of these options consist of the £ = 2, m = +2 modes only, 
with one tapered in time to remove the merger, and the 
other including the full inspiral-merger-ringdown signal. 
The midpoint of the taper corresponds to the time when 
the signal reaches the frequency of the innermost stable 
circular orbit (ISCO) frequency of a test particle orbit- 
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TABLE II: Ratio of the waveform-model results for median variance of all the extrinsic parameters for two sets of comparable- 
mass physical systems. The "Numerator" and "Denominator" columns indicate the models compared in constructing the ratios 
for that row. The models vary by the harmonic content of the waveforms and by whether the merger is included (full) or not 
(ISCO). For the systems considered here, the fractional loss in estimated precision from ignoring the final merger is comparable 
to the corresponding loss in SNR and has a greater impact than ignoring higher harmonics. The significance of the higher 
harmonics is lower when full models are considered, as compared with ISCO-terminated models. The actual median fractional 
variances for all cases are given in Table ULTl 



mg a Schwarzschild black hole, / ISCO = c 3 / (6 3/2 ttGM). 
Much of the previous systematic work on parameter un- 
certainties with LISA observations has applied waveform 
models similar to the (£ = \m\ = 2, ISCO) option. More 
recent work has included higher harmonic content 
|35|. Our other two waveform options include modes up 
to I < 4, where one case is again tapered to remove the 
merger, and the other includes the complete signal. 

The top row for each system included in Table ILTl shows 
the ratio of parameter uncertainties with and without 
the merger when higher harmonics are included. In each 
case the inclusion of the merger increases the SNR by 
roughly a factor of 3. In general terms, if the informa- 
tion contained in the merger waveform is equally rich, 
as compared with the inspiral waveform, the uncertain- 
ties should decrease by a similar factor. This is generally 
the case for most parameters with the mass uncertainty 
showing the least improvement and the sky-position, po- 
larization and coalescence-time uncertainties improving 
most. 

The second and third rows for each system in Table 
ILT1 summarize the improvement in uncertainties result- 
ing from the inclusion of higher-harmonic content in the 
waveforms. The second row shows the improvement in 
uncertainty when the full waveform is included in each 
model, while the third row shows the effect of including 
the higher harmonics with waveform models terminating 
at ISCO. The latter comparison is roughly similar to pre- 
vious considerations of the impact of higher harmonics 
[3lTj33j . Comparing the second and third rows provides 
some indication of the independence of information in 
the higher-harmonics and in the post-ISCO merger. For 
most parameters the marginal effect of including higher 
harmonics is not as great when the full-length waveforms 
are considered as it was for ISCO-terminated waveform 
models. 

In Figs. [3J and SI we show histograms of our results 



for the four waveform model options for the q = 1 and 
q = 1/2 systems, respectively. The histograms in Fig. [3] 
agree qualitatively with those presented in our prior work 
[9|. Quantitatively there is some disagreement, which 
may be attributed to several factors. Chief among these 
are the increased duration of the signal (~10 days in the 
prior work as opposed to ^3 months for a comparable 
mass in this work) and an error in the prior code that 
omitted a factor of the TDI cadence in the parameter 
uncertainty estimates. 

For both the equal-mass case in Fig. |3J and q — 1/2 in 
Fig. 21 we see a clear improvement in the level of measure- 
ment precision one can expect by including the merger 
waveform. It appears that, in particular, the parameter 
t c is localized extremely well in both Figs. [3] and 2] for 
cases that include the merger, relative to the timing pre- 
cision without the merger. Indeed, for every mass ratio 
the inclusion of the merger is estimated to result in un- 
certainties in t c that are an order of magnitude or more 
smaller than the smallest gravitational wave half-period 
reached by the signal waveforms, which is the shortest 
time interval over which the signal will contain informa- 
tion content. This dramatic improvement in timing ac- 
curacy can be heuristically explained by noting that the 
merger provides a sharp feature that can be well local- 
ized. 

The total system mass, M, is essentially insensitive to 
the inclusion of the merger or the presence of higher har- 
monics, and appears to depend entirely on the number 
of inspiral cycles, as was anticipated in [261 ] . For this rea- 
son, we instead show SNR -1 in Fig. 0] and in subsequent 
histograms. This quantity is useful, as it shows the de- 
gree of relative improvement in parameter measurement 
that can be explained by an increase in SNR alone. 

The precision of the luminosity distance Dl and po- 
larization phase ip measurements improve by roughly an 
order-of-magnitude over the quadrupolar inspiral case as 
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cither the merger or higher harmonics are added indi- 
vidually. When both features are added simultaneously, 
the improvement is "only" a factor of ~ 30, suggesting 
that some of the information added by the two features 
is common. 

The inclination t and orbital phase constant <p a show 
qualitatively different behavior depending on what addi- 
tional physics is added to the waveform model. For both 
mass-ratios, we see that the addition of higher harmonics 
dramatically reduces the long tail of large uncertainties 
for the parameters in the worst cases of the quadrupolc- 
only results. The addition of the merger, on the other 
hand, results in an improvement of the most precise de- 
terminations of i and 4> , with less effect on the uncer- 
tainty of the least accurate parameter sets. Unlike the 
results for luminosity distance and polarization phase, 
the overall improvement when both merger and harmon- 
ics are included is closer to the product of the individual 
improvements, indicating that the additional information 
brought by each is independent. 

For the sky angles (the ecliptic latitude j3 and longitude 
A) the phenomenology of the response to including higher 
harmonics is roughly reversed from that seen in the in- 
clination and orbital phase angles. In both Figs. [3] and IH 
inclusion of higher harmonics most significantly improves 
the smallest uncertainties in the distribution, with less 
impact on the largest uncertainties. The addition of the 
merger, however, appears to be more complicated. For 
both mass ratios the £ = \m\ = 2 uncertainty distribu- 
tions for both sky angles appear to uniformly improve. In 
the £ < 4 distributions, the addition of the merger shows 
relatively more improvement in the least-accurate side of 
the distribution. The reduction in uncertainty obtained 
by adding both features to the waveform model shows less 
independence than seen with i and <$><,. For the q = 1/2 
in particular, there is relatively little benefit to adding 
the higher harmonics once the merger has been included. 



B. Systems with different total masses 

With rcdshiftcd mass 1.33 x 10 6 M Q , the frequency 
of the inspiral-mcrgcr transition in the signals we have 
studied so far occurs near the optimal region of LISA's 
sensitivity band. Varying the mass shifts this transition 
frequency (in inverse proportion), thus changing LISA's 
relative sensitivity to the inspiral and merger-ringdown 
signals. 

In Fig. [5j we compare three cases, all with a mass ra- 
tio q = 1/2 at a redshift z = 1, and with total masses of 
M = 1.33 x 1O 5 M , 1.33 x 1O 6 M , and 1.33 x 10 7 M Q , 
chosen in part so that the heaviest case can be compared 
to the results in ||. In this and subsequent figures, we 
do not compare results for the measurement of the mass 
due to the fact that our signal duration is constant in 
M, and therefore the lightest systems do not fully span 
LISA's band. This was essentially due to computational 
constraints on the signal length, which we intend to im- 



prove upon in future work. 

Scaled in units of seconds, we see the best t c estimates 
for the mid-mass case, which merges closest to LISA's 
most sensitive band, and therefore has the largest signal- 
to-noise ratio (SNR). On the other hand, if we were to 
rescale the curves to measure precision against the time- 
scale of the source physics, M, then the largest systems 
would be seen as most precise. 

For the sky position angles the middle mass case out- 
performs the others by a factor of 2-3, with a broad dis- 
tribution for the highest mass case. For all other pa- 
rameters, the lowest mass is easily the worst performer, 
with the mid-mass system marginally outperforming the 
largest mass case. In prior investigations that were 
limited to the inspiral a more precipitous drop in per- 
formance occurs for systems with masses approaching 
1O 7 M . This is simply a result of the absence of signal, 
as for such large masses the portion of the total signal 
that occurs in band for LISA is increasingly dominated by 
the merger, so that no signal is present when the merger 
is excluded. This effect is exacerbated in studies that 
employ a more severe low-frequency cut-off in the LISA 
sensitivity. 



C. Systems with different mass ratios 

We have also examined results for mass ratios other 
than q = 1/2. In Fig. |6l we compare three different 
mass ratios, q = 1/2, q = 1/4, and q = 1/10, with all 
three cases again corresponding to a total system mass 
of 1.33 x 10 6 M o at a redshift z = 1. 

Varying the mass ratio has surprisingly little effect on 
the uncertainties. We see that the inverse SNR shows 
more variation than the parameter uncertainties. This 
would seem to suggest a balance between the importance 
of the total signal power and the fraction of that power 
contained in higher harmonics. To the extent that the rel- 
atively small differences among the three cases for most 
parameters are statistically meaningful, the q = 1/10 
case is the worst performer by a small margin for all 
parameters except t c , with an insignificant difference be- 
tween the q = 1/2 and q = 1/4 cases. 

D. Comparing results 

In Table Hill we summarize our results by quoting the 
median parameter uncertainties for all of our data, as 
well as quoting results for comparable cases from the lit- 
erature. We note that, as we do not include the mass 
ratio in our covariance calculation, we are unable to con- 
vert to uncertainties in the chirp mass, Ai c , and reduced 
mass ratio, which are used in [32j and in most of the 
literature. Furthermore, without explicit knowledge of 
the covariance between these parameters in the available 
publications, we are unable to convert the results in the 
literature into uncertainties in the total mass, so we leave 
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FIG. 3: Uncertainty histograms for q = 1 at redshift 2 = 1, corresponding to a total system mass M = 1.33 x 1O 6 M , calculated 
using a full inspiral-merger-ringdown waveform with harmonics £<4 (solid), an inspiral waveform truncated at the ISCO as 
described in Sec. IIII Al (dashed), a full waveform including only quadrupole (I = 2, m = +2) modes (dash-dotted), and an 
inspiral waveform including only quadrupole modes (dotted). All histogram bins are normalized by the total number of cases 
and are expres sed as percentages. 
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FIG. 4: Uncertainty histograms for q = 1/2 at redshift 2 = 1, corresponding to a total system mass M = 1.33 x 10 6 M Q , 
calculated using a full inspiral-merger-ringdown waveform with harmonics I < 4 (solid) , an inspiral waveform truncated at the 
ISCO as described in Sec. IIII A1 (dashed), a full waveform including only quadrupole [L = 2, m = +2) modes (dash-dotted), 
and an inspiral waveformincluding only quadrupole modes (dotted). All histogram bins are normalized by the total number of 
cases and are expressed as percentages. 
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FIG. 5: Uncertainty histograms for q = 1/2 at z = 1, corresponding to a total system mass M = 1.33 x 10 5 M Q (solid), 
1.33 x 10 6 M Q (dashed), and 1.33 x 1O 7 M (dash-dotted). All histog ram bins are normalized by the total number of cases and 
are expressed as percentages. 
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FIG. 6: Uncertainty histograms for M = 1.33 x 1O 6 M0 at z = 1, for mass ratios q = 1/1 (solid), 1/2 (dashed), and 1/4 
(dash-dotted). All histogram bins are normalized by the total number of cases and are expressed as percentages. 
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the mass out of our comparison in Table IIIII We do in- 
clude a comparison of the sky position, calculated using 

CD- 

Of particular note is the discrepancy between our re- 
sults, and the results found in Q. Specifically, from their 
Fig. 1, their median latitude and longitude uncertain- 
ties are 0.046 deg and 0.057 deg, respectively, or ~ 3 
arcmin as stated in their abstract. However, when run- 
ning as identical a case as possible with ten cycles prior to 
merger of a system with q = 1/2 and M = 1.33 x 10 7 M Q , 
and using the A and E channels only, we arrive at me- 
dian estimates of 0.39 deg and 0.63 deg for the latitude 
and longitude, respectively. This represents an order-of- 
magnitude disagreement. We note that for a total mass 
of M = 1.33 x 1O 6 M0, the median latitude and longi- 
tude uncertainties for all mass ratios was within a factor 
of a few of the localization claimed in Q, with q = 1/2 
providing the best localization with median latitude and 
longitude uncertainties of 0.09 deg (5 arcmin) and 0.18 
deg (11 arcmin), respectively, and with 10% of the cases 
in that ensemble being localized at the ~ 1 arcmin level. 



E. Accumulation of information with time 

In actual observations of black hole binaries with LISA, 
information about the system will be progressively un- 
veiled over time. In particular some estimate of the 
system parameters may be available in advance of the 
merger observation. As the system approaches merger 
the uncertainties of these estimates are expected to de- 
crease sharply [l4| . This real-time development is es- 
pecially important in planning multi-messenger observa- 
tions. How and when sky position estimates improve 
as the coalescence proceeds may impact the instruments 
operational requirements including how frequently data 
downlinks are required and in planning protected observ- 
ing periods near the moment of merger. In turn, these 
operational requirements may influence details of LISA's 
instrumental design. 

In order to compare the evolution in measured param- 
eter precision for different systems, we have calculated 
the parameter uncertainty for waveforms whose ends are 
"turned off" via windowing as described earlier. In this 
section, all the specified times correspond to the time 
at the mid-point of the applied taper. This procedure 
is analogous to a realistic procedure for measuring pa- 
rameters from progressively longer segments of real-time 
data. 

In Fig. we compare the uncertainty in ecliptic lati- 
tude (5 and longitude A for the equal-mass waveform and 
the g=l/2 waveform, both with higher harmonics £<4 
and restricted to quadrupolar modes (£ = 2,m = +2). 
The linear-appearing decrease seen in both panels indi- 
cates that, over the last several hours before merger, our 
estimates for uncertainties in the sky position angles are 
roughly proportional to the time remaining before merger 
until a couple minutes before "merger" (which we have 



defined as the moment at which the \£\ = 2 mode ampli- 
tude peaks). For the cases studied, the dominant ring- 
down radiation period is about 80 s, roughly setting the 
scale at which the linear trend levels off. Note that, in 
some cases the parameter uncertainties may continue to 
improve after "merger" , drawing on information in the 
ringdown radiation. The lower pair of curves in each 
panel are based on the waveform model including the 
higher harmonics. Consistent with the discussion in Sec. 
IIII A[ including the harmonics continues to be valuable 
even late in the observation, after the merger is recorded. 



Fig. |8] shows another comparison of latitude uncer- 
tainty, for the same systems compared in Fig.[SJ Because 
the parameter being varied is the total system mass, and 
the mass rescales time, we compare these results using 
times measured in M and in seconds. Because these 
signals are simply mass-rescalings of the same signal in 
naturalized units, and therefore have identical harmonic 
content relative to their quadrupolar content, the main 
factor for differences is the frequency band spanned by 
the signal, and where that band falls relative to the most 
sensitive band of the detector. The lowest-mass case, 
M = 1.33 x 10 5 M Q , has the largest number of cycles in- 
band, and therefore performs best at early times. By the 
time it merges, however, the signal is chirping at frequen- 
cies much higher than the most sensitive band for LISA, 
so the the contribution after ISCO is negligible for this 
case. The mid-mass case, M = 1.33 x 10 6 M©, is out- 
performed by the lowest-mass case at early times. How- 
ever, because it merges in LISA's most sensitive band, 
the contribution approaching ISCO and running through 
the merger and ringdown is far greater than the other 
cases. Indeed, by the time the full signal has been in- 
cluded, this case yields a more accurate estimate than the 
lowest-mass case by a factor of ~ 2. The largest mass, 
M = 1.33 x 10 7 M Q , has the fewest cycles in band, so it 
yields the lowest precision at early times. However, it too 
has a substantial gain in SNR, relative to the SNR of its 
inspiral, in the late inspiral through the merger and ring- 
down, so it too makes gains in precision relative to the 
lowest-mass case, although unlike the M = 1.33 x 10 6 M Q 
case, it does not fully "catch up" , and remains the worst 
performer of the three. 



We compare latitude uncertainty for four different 
mass ratios in Fig. |H1 q = 1, q = 1/2, q = 1/4, and 
q = 1/10. This comparison again shows a trade-off be- 
tween the number of in-band cycles and the signal power. 
Because radiation reaction is weaker for more disparate 
mass ratios, the q = 1/10 yields the highest precision at 
early times, despite having significantly less power (the 
SNR scales as rj for the inspiral, and as r] 2 for the merger 
[3, 0)- Sky position uncertainty for the smaller-^ cases 
decreases more slowly over most of the last day than the 
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"signal duration is limited to 10 cycles, for comparison to Q 
''estimated from histograms in Fig. 2 of [Til 

C "R" indicates that the amplitude was restricted to the leading 
order term. 

^estimated from the results in Table II of [2^1 

TABLE III: Median fractional variance of all the extrinsic parameters for the cases investigated in this paper, as well as results 
from the literature for comparable systems. All angles are measured in degrees, and time is measured in seconds. We separate 
our results into cases where all available information has been included (top portion), and where some information has been 
suppressed for testing and analysis (bottom portion). For literature results, "A PN" refers to the post-Newtonian order of the 
model used. All studies set at a fixed source distance of z = 1, albeit with slightly different cosmological parameters. The 
results in this work correspond to ~ 10 6 M of observation, (~ 3 months for M = 1.33 x 10 6 Mq) unless otherwise noted, while 
the results from the literature correspond to 1 year of observation. 




FIG. 7: Comparison of the median ecliptic latitude (left panel) and longitude (right panel) uncertainties as a function of time 
before the merger for a mass M — 1.33 x 10 6 Mq at z = 1, for mass ratios of q = 1 and q — 1/2. For each mass ratio, we 
compare the uncertainty for a waveform including harmonics £<4 (q — 1: solid, circles; q = 1/2: dash-dotted, down arrows), 
and for waveforms including only the quadrupole [i = 2, m = +2) modes (q — 1: dashed, squares; q = 1/2: dotted, up arrows). 
The vertical dash-dotted line corresponds to the Schwarzschild ISCO. 
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FIG. 8: Comparison of the latitude uncertainty as a function of time before the merger for a mass ratio q = 1/2 at z — 1, 
with a total system mass M = 1.33 x 1O 5 M (solid, circles), M = 1.33 x 10 6 M o (dashed, squares), and M = 1.33 x 1O 7 M 
(dash-dotted, down arrows). In the right panel, the time axis is scaled in units of M, the physical time-scale of the merger 
process, rather than seconds. The vertical dash-dotted line corresponds to the Schwarzschild ISCO. 



near-linear rate seen for the equal-mass case. By ~ 20 
minutes before merger, approaching ISCO, the median 
uncertainties in j3 arc roughly the same for all mass-ratios 
shown. At late times, the power content becomes a more 
dominant factor in further decreasing the uncertainty in 
/?. The equal-mass case contains more signal power in 
the merger. The q = 1/2 and q = 1/4 cases are nearly 
optimal, retaining some of the merger signal strength but 
perhaps benefiting more from stronger harmonic content 
at late times [l6[. Overall final sky-position error esti- 
mates are nearly flat for 1 > q > 1/4 (see Tabic . 
By q = 1/10 the merger signal power is significantly di- 
minished. While there are still improvements in position 
estimates after merger, they are notably smaller than 
those in the other cases. 
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»q = 1/2 
f q = 1/4 
*q = 1/10 
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FIG. 9: Comparison of the latitude uncertainty as a function 
of time before the merger for a mass M = 1.33 x 10 6 Mq at 
z = 1, and mass ratios of q = 1 (solid, circles), q — 1/2 
(dashed, squares), q = 1/4 (dotted, down arrows), and q — 
1/10 (dash-dotted, up arrows). The vertical dash-dotted line 
corresponds to the Schwarzschild ISCO frequency. 



IV. CONCLUSION 



We have investigated the precision with which black 
hole binary system parameters can be measured from 
LISA observations including merger waveforms in the 
analysis of nonspinning binaries with moderate mass ra- 
tios (<7>1/10). We have further studied how the ex- 
pected performance depends on mass ratio and total sys- 
tem mass, and the impact of including or neglecting the 
merger signal and higher harmonics. The luminosity dis- 
tance and the polarization phase depend on both the in- 
clusion of the merger and the presence of higher har- 
monics, although the improvements from including these 
two elements are not independent. The inclination, the 
orbital phase constant, the ecliptic latitude and ecliptic 
longitude also depend on both the merger and the har- 
monic content. For these parameters the improvements 



resulting from including both the merger and higher har- 
monics arc essentially independent. 

For comparable-mass systems near 10 6 M Q , ignoring 
the merger reduces the SNR by a factor of ~ 3 and re- 
sults in a similar loss of median precision in parameter es- 
timation, even more so for the sky position estimates. For 
sky position, ignoring the merger results in a more sig- 
nificant loss of precision than ignoring higher harmonics. 
Parameter estimates are roughly independent of mass ra- 
tio through 1 > q > 1/4, for sky position in particular, 
though for smaller mass ratios q < 1/10 the precision 
begins to decrease. For q = 1/2, the best parameter esti- 
mates are obtained for systems near 1O 6 M0, which merge 
in the middle of LISA's sensitivity band. Decreasing the 
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mass by an order of magnitude to ~ 1O 5 M0 results in 
a precision loss of roughly a factor of 5 with a dimin- 
ished relative contribution from the merger. Increasing 
the mass to ~ 10 7 M© results in a similar loss. Though 
we have left out the effects of spin, including precession, 
our median sky position precision estimates are similar to 
those obtained with precession, but ignoring the merger 
fl4| . Each method locates the systems in the sky within 
O(10 arcmin). Our best cases (~ top 10%) are local- 
ized at the level of 0(1 arcmin). We estimate that LISA 
will usually be able to locate larger mass systems (near 
10 7 M o ) quite well, in some cases better than systems 
with masses near 10 5 Mq, and far better than earlier es- 
timates based on inspirals alone. Our results for these 
more massive systems do not, however, reproduce the 
preliminary (but widely discussed) extraordinarily pre- 
cise sky localization results found in [||, though we do 
achieve such high precision for the ~ 10 6 M Q systems, 
where both the inspiral and merger are in-band and can 



contribute. For equal-mass systems near 10 6 M Q the sky 
angle estimates improve over the last several hours up to 
a few minutes before merger in rough proportion to the 
time remaining before merger. 
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